Theory for Raman superradiance in atomic gases 
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A mean field theory for Raman superradiance (SR) with recoil is presented, where the typical 
SR signatures are recovered, such as quadratic dependence of the intensity on the number of atoms 
and inverse proportionality of the time scale to the number of atoms. A comparison with recent 
experiments and theories on Rayleigh SR and collective atomic recoil lasing (CARL) are included. 
The role of recoil is shown to be in the decay of atomic coherence and breaking of the symmetry of 
the SR end-fire modes. 



I. INTRODUCTION 

Superradiance, first proposed by Dicke 0], is the en- 
hanced radiation from a collection of coherently decaying 
dipoles. It has been studied extensively theoretically (see 
review and references therein) and has been observed 
in many different systems, including thermal gases 0, 
and Bose-Einstein condensates (BECs) 00 El El- There 
is mostly agreement now on the fact that the collectiv- 
ity is responsible for superradiance, which is the same 
no matter whether the medium consists of Bosons or 
Fermions H, ^ n tne case °f BEC, collectivity can 
be observed as matter wave stimulation, or "Bosonic en- 
hancement" ■ BEC is unique in that there is negligible 
Dopplcr broadening and the recoil momentum is mea- 
sured easily and in fact was recently used to demonstrate 
BEC superradiance 0,01 IS El- In particular, superradi- 
ance can be described by "collective atomic recoil laser" 
(CARL) equations in the bad cavity regime pi flOL fll| . 
Collective gain can be observed with CARL in the sense 
that it depends nonlinearly on the density [T3 | and thus 
does not occur for atomic densities below a certain criti- 
cal value . 

Most experiments on superradiance were done using 
pulsed pump lasers to "instantaneously" invert a two- 
level system. The quantum stage of superradiance, where 
the radiation field builds up from vacuum fluctuations, 
can then be modeled to start only after the pump laser is 
turned off 0. For this case, pump lasers obviously have 
to be strong; at the same time, experiments done with 
BECs use only weak pump fields. We therefore call the 
one strong pump superradiance and the other weak pump 
superradiance. In the latter case, the quantum stage hap- 
pens while the pump field is still on. In this article, we 
will focus on weak pump superradiance. Note that in 
this case the maximum instantaneous superradiance rate 
is limited by the pump laser intensity, while for strong 
pump superradiance no such limitation exists. 

Mostly, earlier research concentrated on so-called 
Rayleigh superradiance E E E which 
happens for transitions between different center of mass 
(cm.) states while the internal state remains un- 
changed [2jj. We will here discuss Raman superradi- 
ance, where there are two different internal ground states 



for the pump and the superradiant transition. Recoil 
and different cm. states are taken into account here as 
well, but are, as we will show, of lesser consequence. It 
turns out that Raman superradiance otherwise follows 
the same basic patterns as Rayleigh superradiance. Al- 
though superradiance with Raman pumping has been 
analyzed in Ref. [2l|, the recoil effect was ignored and 
the Raman pumping time was assumed to be short com- 
pared with the superradiance time. It will be shown in 
this paper that recoil induces the decay of Raman coher- 
ence and may make the superradiant modes asymmet- 
ric. In Ref. |22| an incoherent cw pump laser was con- 
sidered numerically, also leading to superradiance. Re- 
cently, M. M. Cola, et al. j2^| presented a quantum theory 
to describe the Raman superradiance experiments with 
BECs 0, IH 0- In comparison, our analysis can be ap- 
plied to both thermal atoms and BECs with emphasis on 
the effect of recoil. We also discuss the connection with 
CARL using stability analysis. In addition, we consider 
the asymmetry of superradiant modes as the pump laser 
setup is changed which helps to understand the underly- 
ing physics of superradiance. 

This paper is organized as follows: In Sec [n] we derive 
the dynamical equations to describe Raman superradi- 
ance. These equations are used to analyze the stability 
conditions in Sec. IIIII Numerical calculations in compar- 
ison with experiments are included in Sec lIVI Discussion 
and conclusion follow in Sec. 



II. MODEL 



We consider a three-level A-type atomic system with 
excited state |1) and two ground states |2) and |3) 
(Fig. When the detuning of a pump laser is much 
larger than both, its Rabi frequency and the maximum 
Rabi frequency of the superradiant field, the interac- 
tion picture Hamiltonian of this system under dipole and 
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rotating- wave approximation reads |2J, 1251 l26| 

+ E hgl^-^+^s + H.c. 

+ E nglge-^- % ^iap> 3 + H.c. , (1) 



with coupling constants g* ~ 



hkc 
2e V 



£2 • di2yr and 



'3 k = l \J ^oF e " 3 ' ^ 13 ^2~' ™ wnat follows assumed to be 
real, ei is the polarization direction. While the cm. 
Hamiltonian is H cm = — |^V 2 with m being the mass 
and ?i being the Planck constant, the Hamiltonian of the 
optical fields is Hf = J2q? l( l ca q* 2 a q2 + ^ca^a,^, where 
o-q2 {o-q3) is the field annihilation operator for the tran- 
sition between |1) and |2) (|3)), and q the momentum of 
the radiation field. is the atomic field operator. V 
the quantization volume, £1 and A the pump field Rabi 
frequency and detuning, 5 3 the ac-Stark shift due to the 
pump laser, dyi idvi) the dipole moment between |1) and 
1 2) (|1) and |3)). In Eq. Q, the second line describes the 
Rayleigh transition, the third line the Raman transition. 
The ratio of <72,<f/ 1 9s,q determines the branchingratio be- 
tween Rayleigh and Raman superradiance While 
the Ra yleig h superradiance has been studied extensively 
13 H3i this paper will focus on Raman superradi- 
ance. 




FIG. 1: Center of mass manifolds associated with three- 
internal-state atomic system. State |2) is the one-particle 
state of the initial BEC. Both Rayleigh transition and Raman 
transition are present with Raman field a^p and Rayleigh field 
a^. The pump laser Rabi frequency Q is much smaller than 
the detuning |A|. 

Using Fock representation, Eq. can be written as 



H/h 
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-/i.c.)+E< 



i q2 a q 2 
(2) 

where bjk (j = 2, 3) annihilates an atom in state \j) 
with momentum k and energy w k = h 2 k 2 /2m. b 2k = 
b2,k+k - q , u q = cq - uj + oj 23 + S 3 , uj = ck and u 23 is 
the atomic energy difference between |2) and |3). For 
simplicity, the vector arrows from q and fc have been 
dropped here. We assume that g 2 , q ~ 32 only weakly 



depends on q for the relevant range of modes. From this 
form it is obvious that the total population on |2) and 
1 3) is conserved. The matter wave mode b 3k is coupled 
to different b 2k f° r different optical modes q. When the 
detuning is large, however, collective linewidth or multi- 
ple scattering can be neglected and we can drop the 
coupling between different modes. In this article, we also 
neglect the depletion of BEC due to other modes. Raman 
transitions in different directions can thus be considered 
independently. From Eq. J5J), Maxwcll-Bloch equations 
can be derived: 



dt 

d 

Jt P ' hk 



-iuj k A - ig 2 E Pkk - kA, 



(3a) 



-%{u k -u li ) Pik -ig2(l-2N k )A, (3b) 
i(92P kk A- c.c.) , (3c) 



where A = (a q2 ), p kk = (b 2k h k ), N k = (b 3k b3k), with N k 
being the number of atoms in state |3/c). n is the effec- 
tive radiation field decay rate, if we neglect propagation 
in the mean field approximation p?! l29j . This approx- 
imation works well when the medium is optically thin 
at the pump frequency, which is the case here since the 
pump field is far detuned from resonance. With L and D 
the length and diameter of the medium and A the wave- 
length of the superradiant transition, the Frcsnel number 
F = jj\ gives approximately the number of modes that 
fit in the medium in axial direction. If it is around or 
bigger than 1 as in the experiments then k = c/2L 

for axial modes (also called "end- fire modes" 0), which 
are the modes having largest gain for superradiance, and 
K off > <ni(j?+l) f° r off-axial modes J29^ . It will be shown 
in Sec. 11111 that in experiments 0, 0], k dominates over 
all the other relevant characteristic rates and therefore 
makes the end-fire modes most likely to superradiate. In 
the following, we assume all superradiant modes to be 
axial. 

In a BEC, only the k = state is present, and thus 
Eqs. © become 



dt 
d 

Jt Pm 
d 

dt N ° 



-iu k A - ig 2 Np 0o - kA, (4a) 
iu rPm -ig 2 {l-2N a )A (4b) 
i{92Pm A ~ c.c.) (4c) 



where p 5 o = (^o^ 3 - ) witl1 & 2,o = &2,fc -g) the recoil en- 
ergy Tiuj r = h 2 (ko — q) 2 /2m, and N the total number of 
atoms in the system. Because we assume k to be very 
large it follows from Eq. (|4a|l that 



(5) 



-ig 2 Np m /K 
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Substituting A into Eqs. i|4bl4c[) . we arrive at 



dr 
d 



Poo 



No 



g 2 (2N - l)poo 

2|„ |2 



IVrPdO - 
-2#2 I Poo 



(6) 



Here, we scale the time such that r = Nt. It is therefore 
obvious that the timing of the resulting process scales 
with 1/iV, in the same way as in traditional superradi- 
ance Q- From Eq. © we know that the output field 
amplitude A is proportional to N and thus the intensity 
is proportional to N 2 . These are typical characteristics of 
supcrradiance. Note that without recoil io r = 0, Eqs. © 
are completely equivalent to Eqs. (6.36) of @], which de- 
scribe standard superradiance: a radiation cascade down 
the pseudo-spin ladder from J z = N/2 to J z = —N/2, 
giving a hyperbolic secant solution for the dependence 
of the upper-level population on time Q. For BEC, the 
term with u) r , which is due to recoil, only contributes to 
the phase evolution of Raman coherence, not to its decay, 
while for thermal atoms recoil docs induce the decay of 
Raman coherence, as discussed in the next paragraph. 

For thermal atoms, Eq. (|3b|) describes quantum diffu- 
sion as well as generation of Raman coherence . In par- 
ticular, the term —i(ui k — u k )p kk in Eq. i|3bfl shows that 
coherence stored in different levels experiences quantum 
diffusion, since the term will have different values for 
different k. To understand how the quantum diffusion 
works, we assume Raman coherence has been generated 
uniformly for all levels, which means p kk (0) = p(0)p k 
with p(0) being the coherence for one level and pt the 
probability distribution of atom at level fc. If we set 
field amplitude A to zero, the solution of Eq. I|3b(l is 
Pkk(t) = p(0)pfc e _I ^ fc_wt '*. The coherence is then 
P(t) = EkPkkit) = pWEkPke- 1 ^-^- To pro- 
ceed, we need to specify p k at temperature T. Here, 
either Bose-Einstein distribution for Bosons or Fermi- 
Dirac distribution for Fermions are appropriate. For 
simplicity, however, we assume Lorentzian distribution 
Pk = ^ k 2+ s 2 with Sp 2 /2m = kgT/2, which describes 
the atoms well even at sub- recoil temperatures |3(3] • The 
summation can be approximated by an integral and it 
follows that the Raman coherence decays exponentially 



Pit) 



5p 



7T k 



Sp' 



■p- kk (t)dk = p(0)c 



(I) 



where T = 2k ^k B T/ m sin 6</2, lot = 2hk 2 sin 2 9/2 and 

6 is the angle between q and kg. It is clear now that 
the decay rate T depends on the pump laser direction 
fco relative to the superradiant pulse direction q. Thus, 
Eq. I)3bj) can be rewritten as 



— p = (iuj r - T)p - ig 2 (l - 2p 33 )A 
dt 



(8) 



where P33 = ^ fe N k is the population in state |3). If we 
would use a Gaussian rather than Lorentzian density of 



states, the inverse 1/e decay time would be ^/2T rather 
than r. As an example, for Rb at the Doppler limit 
temperature of 143^, V = 1.35 x 10 6 a _1 . 

Comparing the result for thermal atoms in Eq. |JSj with 
the result for a BEC in Eq. Ij4b() . we see that thermal 
distribution contributes additional coherence decay, oth- 
erwise these equations are the same as expected. We can 
therefore generalize the results to 



dt 



A = -iuo k A - ig 2 Np - kA, (9a) 
iu) T p - ig 2 {l - 2p 33 )A - n R p (9b) 
P33 = i(92pA-c.c.) , (9c) 



dt 
d 

~di' 



where the total coherence decay kr = k' r + T. k! r can 
be introduced phenomenologically to contain collisions, 
magnetic gradients, etc., and T = for BECs. These 
equations are now analogous to Eqs. (13-15) in Ref. |23j . 
but can be applied to both, BEC and thermal atoms. 



III. LINEAR STABILITY ANALYSIS 

In this section, we will determine the necessary condi- 
tions for Raman supcrradiance to ha ppen, w hich is easi- 
est using linear stability analysis j^Fj, |3l], H3] • Obviously, 
A = 0, p = 0, and P33 = 1 give a stationary solution of 
Eqs. ©. Rewriting Eqs. © for A = + 8 A, p = + 8p, 
and P33 = 1 + 8p 33 leads to a two-dimensional linear sys- 
tem with the characteristic equation 

S 2 + (i (uj k - u) r ) + k + kr) S + {—iu>r + kr) k 

-Ngl + u k uj r + iuicKR. = . (10) 

(The third equation is equivalent to zero in this case and 
can be dropped.) In comparison with the cubic insta- 
bility equation for Rayleigh superradiance |20j |. this is 
a quadratic equation. The physical reason for such a 
change is that for a Rayleigh transition, the initial and 
final internal states are the same and thus only atoms 
with different cm. states may contribute to the gain (see 
Eq. (49) of Ref. [2(|); for a Raman transition, the initial 
and final internal states are different and thus all atoms 
contribute to the gain regardless of the cm. states. 

The above quadratic equation has two roots for 5, 5+ 
and SL. Since S is the exponent of the state vector 
[{8A,8p) = (cL4(0),<5p(0))exp St], the zero solution be- 
comes unstable if at least one of S+ or S- has a positive 
real part. The larger real part (let's call the respective 
Si 



root S + = S' + + iSD is therefore defined as instability 
factor. If S' + > 0, then the system is dynamically un- 
stable, from which the threshold pump intensity can be 
derived. From Eq. qiUfl. it can be easily seen that S' + de- 
pends nonlinearly on the number of atoms N. Note that 
nonlinear dependence on TV is the essence of collective 
instability [20$. In the bad cavity regime as is the case 
for the experiments 0,O|j K is large, the system therefore 
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depends linearly on the atomic density and therefore may 
display superradiant behavior. In a good cavity, however, 
K is much smaller and thus the collective gain depends 
on the density nonlinear ly [Til IT^ | . 

Since ur and lu^ can be shown to have only a mi- 
nor effect on S' + under the experimental conditions of 
Refs. 0, yj, we set u>r = lo^ = 0. In this case, the insta- 
bility factor is 

, -(k + kr) + ^(k + k r ) 2 + 4(7Vg| - kk r ) , , 
b +~ 2 1 J 

In particular, for vanishing coherence decay kr = 0, 
the system is unstable and therefore superradiant for 
any pump laser power. This is different from the case 
of Rayleigh superradiance which always has a non-zero 
threshold pump laser intensity [2(j ■ In the case of ther- 
mal atoms, however, kr can be considerable, and the 
threshold pump intensity is quite high in the bad cavity 
limit. This explains why collective gain was not observed 
in Refs. 0, 0|. It should be possible experimentally to 
minimize the decay due to quantum diffusion if the pump 
laser is collinear with the sample. Raman superradiance 
or collective gain might perhaps be observed in this case 
even in thermal atoms. 

When k is much lar ger t han any other frequency in 
Eq. (fit]] ), i.e., k \J Ng 2 , Wk, w r , as in the experi- 
ments 0, 0, the instability factor can be simplified to 

S'+*- {Ngt - krk) (12) 

K 

In this case, S' + is linear in N, which means experiments 
in Ref. [4j,|5( would be purely in the superradiant regime. 
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FIG. 2: Effect of number of atoms and detuning on the evo- 
lution of (a) the intensity \A^ and (b) the population psj, as 
a function of time. Parameters used in the calculations are 
from Ref. Q. The 1/N dependence of the delay time and the 
iV 2 dependence of the maximum intensity can be clearly seen. 
The finite population left in state |3) is due to the decay of 
Raman coherence. 



IV. NUMERICAL SIMULATIONS 

To compare our theory with experiments, we solve 
Eqs. (0 for both BEC and thermal atoms. In the simula- 
tions, we use the initial value o f p( 0) = (2/7V) 1 / 2 , which is 
determined by quantum noise |29t l33l| . Other parameters 
are calculated using the data in Ref. @: 52 = 0.5 x 10 6 
s _1 , k = 1.76 x 10 12 s _1 . oj r is negligible in this context. 

In Fig. [3 we show that the intensity of superradiance 
is proportional to N 2 and the superradiance delay time 
is proportional to 1/N at least as long as there is no 
Raman coherence decay, i.e., if we assume a BEC. The 
numerical delay time 75 — 150[j,S also reproduces well the 
experimental data 0, . 

Figure [3] shows that when the Raman coherence decay 
rate kr is increased, the radiation intensity decreases and 
the delay time increases. This is similar to two-level su- 
perradiance: dipole-dipole interaction decreases the co- 
herence between atoms and thus competes with superra- 
diance. Because of the effective population mixing caused 
by Raman coherence decay there is always a finite num- 
ber of atoms in the Rayleigh lower state |3) at any time 



for a finite kr. For thermal atoms at Doppler cooling 
limit T = U3nK, T = 1.35 x lOS" 1 and the instabil- 
ity factor S' + is smaller than zero and no superradiance 
happens. 

The roles of photon and atomic coherence are inter- 
twined for superradiance. Collectivity can be attributed 
to either photons or atoms, or both. In the case of weak 
pump superradiance, the pulse exits the medium and 
thus decays much faster than the (atomic) Raman co- 
herence. Thus the intensity of the superradiant pulses is 
small, and stimulation of photons by photons is not crit- 
ical in this case. For example, if in the calculation the 
pump laser is turned off before all the atoms have radi- 
ated and then turned on again, superradiance continues 
nearly at the same point it was interrupted. This is true 
for an interruption that lasts longer than the photon co- 
herence time (which is here just the escape time of the 
photons of about 1 ps), but shorter than the Raman co- 
herence time, which is between infinity and 1 ms in our 
simulations. The conclusion is that atomic coherence is 
more important than stimulated emission in this case for 
superradiance to happen. 

It was claimed in Ref. that the output photon num- 
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FIG. 3: Effect of Raman coherence decay rate on the 
evolution of (a) the intensity \A\ 2 and (b) the population P33 
as a function of time, kr > is responsible for a longer 
superradiance delay time and a lower maximum intensity. For 
this figure, N = 2 x 10 6 . Other parameters are the same as 
in Fig. H 



random dots in the simulation of Rcf. [jj show the effect 
of the fluctuations in this case. In general, many modes 
might fire simultaneously as long as the population in 
state 1 3) is not depleted. If the delay time of one mode is 
shorter than the sum of the delay and superradiance time 
of the axial mode, then this mode also supcrradiates. The 
same is obviously is also true for the competition between 
Raman and Rayleigh superradiance 0. 

Let us consider the symmetry of the superradiating 
modes. The two axial modes in opposite directions in the 
experiments 0, 0, [B| show identical behavior: the recoil 
pattern is symmetric. One of the reasons for this is that 
the field decay rate for these modes is the same. Also the 
recoil induced decay Kr is zero for a BEC. This is also 
true for an initially fully inverted two level system 0. 
However, if the pump laser is parallel to the sample axis 
(longitudinal pumping), the recoil induced decay for ther- 
mal atoms can be cancelled if the superradiance mode is 
parallel to the pump laser. This breaks the symmetry of 
the two axial modes and privileges the parallel mode over 
the antiparallel one. In Fig.[3|we see that the mode with 
small kr is stronger than other modes and may suppress 
superradiance for them by depleting the population P33. 
The broken symmetry indicates that the equivalence of 
a three-level system with a far detuned pump laser and 
a two-level system does not hold in this case. Note that 
if the atoms are not fully inverted, the s ymm etry could 
also be broken due to stored coherence j36j. However, 
for BECs, since recoil does not contribute to the decay of 
Raman coherence significantly, two superradiant modes 
would still fire even with longitudinal pump. 



ber, N p , enhances the superradiance N r oc (N r + N p + 1), 
where N r is the number of atoms having superradi- 
ated. Indeed, it was assumed that N p = N r Q, then 
N r oc (2N r + 1). However, no cavity was used in 
which means the average N p is small within the sample 
and can be neglected, as is done in Ref. 0, Note 
that the collecting of photons in the (ring) cavity modi- 
fies the rate [3 Hill ■ But in a high Q cavity, the coupling 
between atoms and field is strong and a perturbation ap- 
proach of Fermi's golden rule as used in Ref. p| may not 
apply. Detailed analysis of this is beyond this paper. 

Now we consider the effect of the field decay rate n on 
determining the direction of the superradiant field modes. 
If the Fresnel number F is around 1 as in the experiments 
0,0 , the decay rate of the off-axial modes is much bigger 
than that of the axial mode, and thus only the axial mode 
supcrradiates. Since the decay rates of the Raman coher- 
ence for the axial modes and their directly neighboring 
modes are almost equal to each other, the field decay rate 
determines the radiation direction. On the other hand, if 
F is much bigger than 1, the off-axial modes do not have 
a decay rate much different from the axial ones, and thus 
they may also superradiatc. In this case, the quantum 
fluctuation stage determines which modes are fired. The 



V. DISCUSSION AND CONCLUSION 

Rayleigh superradiance does not happen without re- 
coil. In comparison to this, recoil is not critical for Ra- 
man superradiance to happen, which means that atomic 
bunching and density grating pictures do not apply for 
explaining Raman superradiance, as they do for Rayleigh 
superradiance. Interference between pump laser and su- 
perradiance output |27j equally does not apply in a case 
where both transitions radiate light with different po- 
larization. We therefore believe that collective effects, 
which might be called Bosonic stimulation in the case of 
Bosons, are the main players in Raman superradiance. 

Interesting is the relationship between Rayleigh and 
Raman superradiance. States with different momentum 
may be considered to be orthogonal |l4j in the same way 
as different internal levels, and thus the Rayleigh transi- 
tion can be looked upon as a Raman transition between 
different motional states [32, l37l l38l |. Indeed, the gain 
coefficients have a similar functional dependence on the 
atomic density [l3[ . In particular, in the case of thermal 
atoms with a pump laser not parallel to the sample axis, 
i.e., with large kr, it can be shown from Eq. that the 
instability factor S' + dep ends linearly on N. This was the 
regime discussed in [32j|32l|3g] m which the Raman tran- 
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sition is considered to be in the (linear) single-atom gain 
regime [39j . Although atom statistics are not critical for 
superradiance [2^, the Fermi momentum kp in Rayleigh 
scattering is replaced by the relative momentum differ- 
ence in Raman scattering, thus the problem with a very 
short coherence time in the case of fermions due to recoil 
might be overcome ||. As is done for Rayleigh superra- 
diance |40j . also atom-atom interaction can be included, 
and will be presented in a forthcoming publication. 

Finally, we would like to differentiate two concepts: 
collectivity and collective gain or collective instability. 
Collectivity means that all atoms in the system con- 
tribute to the same mode ^jji while collective gain or 
collective instability ^3] means that the gain depends on 
the number of atoms N nonlinearly. While the experi- 
ments are in the non-collective gain regime, collectivity 
still plays a major role in Raman superradiance. Raman 
superradiance therefore shows that it is the collective ef- 
fect rather than "Bosonic stimulation" that is responsi- 
ble for superradiance 0, @ • It was claimed that 
if the pump laser makes the two-photon detuning for su- 
perradiant mode zero, and thus the Rayleigh transition 
corresponds to a Raman transition between different cm. 
states there would be a single-atom gain instead of col- 



lective gain |39j. However, we tried to show that even in 
a pure Raman transition, collective gain is still possible 
if a cavity is included. 

To conclude, we developed a mean field theory for Ra- 
man superradiance. Raman superradiance does not nec- 
essarily have an intrinsic threshold for pump laser inten- 
sity even if the decay of the optical field is included. We 
found that recoil induced decay of Raman coherence may 
break the symmetry of the two axial modes if the atoms 
are pumped longitudinally, in which case it is possible 
to realize Raman superradiance even in thermal atoms 
while at the same time it might not be possible to re- 
alize Rayleigh superradiance. We also note that both 
the Rayleigh and Raman superradiance experiments were 
done in the regime where the pump laser is far detuned, 
such as not to populate the excited state. What happens 
in the case of a resonant pump laser is under investigation 
presently. 
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